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We report an improved low-energy extrapolation of the cross section for the process 7 Be(p, 7 ) 8 B, 
which determines the 8 B neutrino flux from the Sun. Our extrapolant is derived from Halo Effective 
Field Theory (EFT) at next-to-leading order. We apply Bayesian methods to determine the EFT 
parameters and the low-energy 5-factor, using measured cross sections and scattering lengths as 
inputs. Asymptotic normalization coefficients of 8 B are tightly constrained by existing radiative 
capture data, and contributions to the cross section beyond external direct capture are detected 
in the data at E < 0.5 MeV. Most importantly, the 5-factor at zero energy is constrained to 
be 5(0) = 21.3 ± 0.7 eV b, which is an uncertainty smaller by a factor of two than previously 
recommended. That recommendation was based on the full range for 5(0) obtained among a discrete 
set of models judged to be reasonable. In contrast, Halo EFT subsumes all models into a controlled 
low-energy approximant, where they are characterized by nine parameters at next-to-leading order. 
These are fit to data, and marginalized over via Monte Carlo integration to produce the improved 
prediction for S{E). 

PACS numbers: 25.20.-x, 25.40.Lw, ll.10.Ef, 21.10.Jx, 21.60.De 


Introduction — A persistent challenge in modeling the 
Sun and other stars is the need for nuclear cross sections 
at very low energies mm- Recent years have seen a 
few measurements at or near the crucial “Gamow peak” 
energy range for the Sun mm, but cross sections at these 
energies are so small that data almost always lie at higher 
energies, where experimental count rates are larger. The 
bulk of the data must be extrapolated to the energies of 
stellar interiors using nuclear reaction models. 

The models available for extrapolation also have limi¬ 
tations. Qualitatively correct models of nonresonant ra¬ 
diative capture reactions, with reacting nuclei treated as 
interacting particles, have been available since the mid- 
1960s [4j . However, these models suffer from weak in¬ 
put constraints and dependence on ad hoc assumptions 
like the shapes of potentials. Developing models with re¬ 
alistically interacting nucleons as their fundamental de¬ 
grees of freedom is currently a priority for the theoretical 
community, but progress is slow, and models remain in¬ 
complete A b initio calculations employing modern 

nuclear forces may yield tight constraints in the future. 

For the ' B e{p, y) 8 B reaction - which determines the 
detected flux of S B decay neutrinos from the Sun - the 
precision of the astrophysical 5-factor at solar energies 
(~ 20 keV) is limited by extrapolation from laboratory 
energies of typically 0.1-0.5 MeV. A recent evaluation |T| 
found the low-energy limit 5(0) = 20.8 ± 0.7 ± 1.4 eV b, 
with the first error reflecting the uncertainties of the mea¬ 
surements. The second accounts for uncertainties in ex¬ 
trapolating those data. It was chosen to cover the full 


variation among a few extrapolation models thought to 
be plausible. Since the differences among S{E) shapes for 
different models were neither well-understood nor repre¬ 
sented by continuous parameters, no goodness-of-fit test 
was used for model selection. 

Halo EFT ffHTS]. provides a simple, transparent, and 
systematic way to organize the reaction theory needed 
for the low-energy extrapolation. The ' Be + p system 
is modeled as two interacting particles and described by 
a Lagrangian expanded in powers of their relative mo¬ 
mentum, which is small compared with other momen¬ 
tum scales in the problem. The point-Coulomb part of 
the interaction can be treated exactly, and the form of 
the strong interaction is fully determined by the order at 
which the Lagrangian is truncated [T6l - iT9l . The coupling 
constants of the Lagrangian are determined by match¬ 
ing to experiment. This is similar in spirit and in many 
quantitative details to traditional potential model or R- 
matrix approaches. However, it avoids some arbitrary 
choices (like Woods-Saxon shapes or matching radii) of 
these models, is organized explicitly as a low-momentum 
power series, and allows quantitative estimates of the er¬ 
ror arising from model truncation. 

The low-energy 5-factor for 'Be(p, q) 8 B consists en¬ 
tirely of electric-dipole {El) capture from s- and d- 
wave initial states to p -wave final states (which dominate 
'Be + p configurations within 8 B). All models are dom¬ 
inated by “external direct capture,” the part of the El 
matrix element arising in the tails of the wave function 
(out to 100 fin and beyond) [U fSO]- Models differ in how 
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they combine the tails of the final state with phase shift 
information and in how they model the non-negligible 
contribution from short-range, non-asymptotic regions of 
the wave functions. 

Halo EFT includes these mechanisms, and can describe 
5(E) over the low-energy region (LER) at E < 0.5 MeV. 
Beyond 0.5 MeV, higher-order terms could be important, 
and resonances unrelated to the 5-factor in the Gamow 
peak appear. Compared with a potential model, the 
EFT has about twice as many adjusted parameters, too 
many to determine uniquely with existing data. How¬ 
ever, calculations of the solar neutrino flux do not re¬ 
quire that all parameters be known: it is enough to de¬ 
termine 5(18 ±6 keV). We fit the amplitudes of recently 
computed next-to-leading-order (NLO) terms [IS] in El 
capture to the experimental 5(E) data in the LER. We 
then use Bayesian methods to propagate the (theory and 
experimental) uncertainties and obtain a rather precise 
result for 5(20 keV). 

EFT at NLO — The EFT amplitude for El capture is 
organized in an expansion in the ratio of low-momentum 
and high-momentum scales, k/A. A is set by the 'Be 
binding energy relative to the 3 He + 4 He threshold, 1.59 
MeV, so A « 70 MeV, corresponding to a co-ordinate 
space cutoff of ss 3 fm. Physics at distances shorter 
than this is subsumed into contact operators in the La- 
grangian. The 8 B ground state, which is 0.1364(10) MeV 
below the ' Be-p scattering continuum HU |22], is a shal¬ 
low p-wave bound state in our EFT: it is bound by con¬ 
tact operators but the wave-function tail should be ac¬ 


curately represented. To ensure this we also include the 
J 71 ' = \ bound excited state of 7 Be in the theory. 7 Be* 
is 0.4291 MeV above the ground state; the configura¬ 
tion containing it and the proton is significant in the 8 B 
ground state HZ]. The large (~ 10 fm) 'Be-p scatter¬ 
ing lengths play a key role in the low-energy dynamics; 
s-wave rescattering in the incoming channels must be ac¬ 
curately described. This also requires that the Coulomb 
potential be iterated to all orders when computing the 
scattering and bound state wave functions m\m- in¬ 
deed Z7 Be Z p a em M p ss kc = 27 MeV while the binding 
momentum of ®B is 15 MeV, so these low-momentum 
scales are well separated from A. We generically denote 
them by k, and anticipate that k/A « 0.2. Since the 
EFT incorporates all dynamics at momentum scales < A 
its radius of convergence is larger than other efforts at 
systematic expansions of this 5-factor P5H5T] . 

The leading-order (LO) amplitude includes only ex¬ 
ternal direct capture. As the 7 Be ground state is | 
there are two possible total spin channels, denoted here 
by s = 1,2. They correspond, respectively, to 3 5i and 
5 52 components in the incoming scattering state, and 
3 P 2 and 5 P 2 configurations in 8 B. The parameters that 
appear at LO are the two asymptotic normalization co¬ 
efficients (ANCs), C s , for the 'Be-p configuration in S B 
in each of the spin channels, together with the corre¬ 
sponding s-wave scattering lengths, a s uansiGE]. The 
NLO result for 5(E), full details of which will be given 
elsewhere [18], can be written as: 


5(E) 


/(E) ^C 2 |5 EC (E;^(E)) + E s 5 sd 


(E; S s (E)) + e s S C x (E- <S a (E)) | 2 + \V ec {E)\ 2 


(1) 


Here, /(E) is an overall normalization composed of final-state phase space over incoming flux ratio, dipole radiation 
coupling strength, and a factor related to Coulomb-barrier penetration HU. 5ec is proportional to the spin-s El 
HU !33 GU external direct-capture matrix element between continuum 'Be-p s-wave and 8 B ground-state wave 
functions. 5cx is the contribution from capture with core excitation, i.e. into the 7 Be*-p component of the ground 
state. Its strength is parameterized by e s . Since 7 Be* is spin-half this component only occurs for s = 1, so £2 = 0. 
Because the inelasticity in ' Be-p s-wave scattering is small HUM] it is an NLO effect. 

Short-distance contributions, 5 sd, are also NLO. They originate from NLO contact terms in the EFT Lagrangian 
[Tg] and account for corrections to the LO result arising from the El transition at distances < 3 fm. The size of 
these is set by the parameters E s , which must be fit to data. 5eC; 5sd, and 5cx are each functions of energy, E, but 
initial-state interactions mean they also depend on the s-wave phase shifts 6 S . At NLO we parametrize 6 S (E) by the 
Coulomb-modified effective-range expansion up to second order in fc 2 , i.e., we include the term proportional to r s k 2 , 
with r s the effective range (see supplemental material) [13 S3 ESI- Finally, T >eg is the El matrix element between 
the d-wave scattering state and the 8 B bound-state wave function. It is not affected by initial-state interactions up to 
NLO, and hence is the same for s = 1,2 channels and introduces no new parameters. This leaves us with 9 parameters 
in all: C 2 2 , ai >2 at LO and five more at NLO: r lt 2 , Ei, 2 , and £i JBj. 


Data — The 42 data points in our analysis come from for the direct-capture 5-factor in the LER: Junghans 

all modern experiments with more than one data point et.al, (experiments “BE1” and “BE3”) [37], Filippone 
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et.ai, [32], Baby et.ai., HE SO], and Hammache et.ai., 
(two measurements published in 1998 and 2001 ) [TO 142] . 
Ref. [T] summarizes these experiments, and the common¬ 
mode errors (CMEs) we assign are given in the supple¬ 
mental material. All data are for energies above 0.1 MeV. 
We subtracted the Ml contribution of the 8 B 1 + reso¬ 
nance from the data using the resonance parameters of 
Ref. [38]. This has negligible impact for E < 0.5 MeV, 
so we retain only points in this region, thus eliminating 
the resonance’s effects. 

Bayesian analysis — To extrapolate S(E) we must 
use these data to constrain the EFT parameters. We 
compute the posterior probability distribution function 
(PDF) of the parameter vector g given data, D , our 
theory, T, and prior information, I. To account for 
the common-mode errors in the data we introduce data- 
normalization corrections, We then employ Bayes’ 
theorem to write the desired PDF as: 


w(9,{£i}\D',T;I) = pr(D|g,{&};T;/)pr(g,{&}|/), 

( 2 ) 


with the first factor proportional to the likelihood: 


lnpr (D\g, {£;};T; I) = c- 


^ [(1 

j=i 




where S(g;Ej) is the NLO EFT S'-factor at the energy 
Ej of the jth data point Dj , whose statistical uncer¬ 
tainty is (jj. The constant c ensures pr ( g , {^}|D; T; I) is 
normalized. Since the CME affects all data from a par¬ 
ticular experiment in a correlated way there are only five 
parameters one for each experiment. 

In Eq. ([2]) pr ( g , {£;}, 1 1) is the prior for g and {£,:}■ We 
choose independent Gaussian priors for each data set’s 
all centered at 0 and with width equal to the assigned 
CMEs. We also choose Gaussian priors for the s-wave 
scattering lengths (ai, 02 ), with centers at the experi¬ 
mental values of Ref. [33], (25, —7) fm, and widths equal 
to their errors, (9, 3) fm. All the other EFT parame¬ 
ters are assigned flat priors over ranges that correspond 
to, or exceed, natural values: 0.001 < Cf 2 < lfm” 1 , 
0 < 7 - 1,2 < 10fm @31145], -1 < ei < 1, -10 < L\ 2 < 
10 fm. We do, though, restrict the parameter space by 
the requirement that there is no s-wave resonance in ' Be- 
p scattering below 0.6 MeV. 

To determine pr (g, {^}|D; T;/), we use a Markov 
chain Monte Carlo algorithm [46] with Metropolis- 
Hastings sampling [47], generating 2 x 10 4 uncorre¬ 
lated samples in the 14-dimensional (14d) g ® {^} 
space. Making histograms, e.g., over two parame¬ 
ters g 1 and $ 2 , produces the marginalized distribution, 
in that case: pr (gi, g 2 \D;T] I) = J pr (< 7 , {^}\D;T; I) 
d£ 1 ... d^dgs ... dgg. Similarly, to compute the PDF of a 
quantity F(g), e.g., S(E;g), we construct pr ( F\D;T;I ) 
= f pr(g,{&}\D;T;I) S(F-F(g))d^ 1 ... d^ 5 dg, and his- 
togramming again sufhces. 



FIG. 1: (Color online.) 2d distribution for Ci (x-axis) and Cf 
(y-axis). Shading represents the 68% and 95% regions. The 
small circle and ellipse are the lcr contours of an ab initio cal¬ 
culation [48] and empirical results [49) . with their best values 
marked as red squares. The inset is the histogram and the cor¬ 
responding smoothed Id PDF of the quantity [C\ +G|] x fm; 
the larger and smaller error bars show the empirical and ab 
initio values. 


Constraints on parameters and the S-factor — The 
tightest parameter constraint we find is on the sum 
Ci~\-C '2 = 0.564(23) fm -1 , which sets the overall scale of 
S(E) [53] . Fig. [l] shows contours of 68 % and 95% prob¬ 
ability for the 2d joint PDF of the ANCs. Neither ANC 
is strongly constrained by itself, but they are strongly 
anticorrelated; the Id PDF of Cf + Cf is shown in the 
inset. The ellipses in Fig. ©show ANCs from an ab initio 
variational Monte Carlo calculation (the smaller ellipse) 
1481 ; 53] and inferred from transfer reactions by Tabacaru 
et al. (larger ellipse) [49j . These are also shown as error 
bars in the inset. The ab initio ANCs shown compare 
well with the present results. (The ab initio ANCs of 
Ref. H sum to 0.509 fm -1 and appear to be in moder¬ 
ate conflict.) Tabacaru et al. recognized that their result 
was ler to 2 a below existing analyses of S'-factor data; a 
1.8a conflict remains in our analysis. We suggest that for 
S B the combination of simpler reaction mechanism, fewer 
assumptions, and more precise cross sections makes the 
capture reaction a better probe of ANCs than transfer 
reactions. 

Fig. [ 2 ] depicts the 2d distribution of Li and e±. There 
is a positive correlation: in S(E) below the 7 Be-p inelas¬ 
tic threshold, the effect of core excitation, here param¬ 
eterized by ei, can be traded against the short-distance 
contribution to the spin-1 El matrix element. The inset 
shows the Id distribution of the quantity 0.33Li/fm —ei, 
for which there is a slight signal of a non-zero value. In 
contrast, the data prefers a positive Z 2 : its Id pdf yields 
a 68 % interval —0.58 fm < L 2 < 7.94 fm. 

We now compute the PDF of S at many energies, and 
extract each median value (the thin solid blue line in 
Fig.©, and 68 % interval (shaded region in Fig.©. The 
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FIG. 2: (Color online.) 2d distribution for ei (x-axis) and 
L\ (y-axis). The shaded area is the 68% region. The inset 
is the histogram and corresponding smoothed Id PDF of the 
quantity 0.33 L i/fm — ei. 
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FIG. 3: (Color online.) The right panel shows the NLO S- 
factor at different energies, including the median values (solid 
blue curve). Shading indicates the 68% interval. The dashed 
line is the LO result. The data used for parameter deter¬ 
mination are shown, but have not been rescaled in accord 
with our fitted {£i}. They are: Junghans et.al, BE1 and 
BE3 22 (filled black circle and filled grey circle), Filippone 
et.al, 25] (open circle), Baby et.al, •.'iff. 1TT| (filled purple di¬ 
amond), and Hammache et.al, [411142] (filled red box). The 
left panel shows Id PDFs for 5(0) (blue line and histogram) 
and 5(20 keV) (red-dashed line). 



PDFs for 5 at E = 0 and 20 keV are singled out and 
shown on the left of the figure: the blue line and his¬ 
togram are for E = 0 and the red-dashed line is for 
E = 20 keV. We found choices of the EFT-parameter 
vector g (given in the supplemental material) that cor¬ 
respond to natural coefficients, produce curves close to 
the median S(E) curve of Fig. [3j and have large values 
of the posterior probability. 

5(20 keV) and the thermal reaction rate — Table [T] com¬ 
piles median values and 68% intervals for the 5-factor 
and its first two derivatives, S'/S and 5"/5, at E = 0 



S (eV b) 

S'/S (MeV -1 ) 

5"/5 (MeV -2 ) 

Median 

21.33 [20.67] 

-1.82 [-1.34] 

31.96 [22.30] 

+CT 

0.66 [0.60] 

0.12 [0.12] 

0.33 [0.34] 

—a 

0.69 [0.63] 

0.12 [0.12] 

0.37 [0.38] 


TABLE I: The median values of 5, S'/S, and S''/S at E = 0 
keV [E = 20 keV], as well as the upper and lower limits 
of the (asymmetric) 68% interval. The sampling errors are 
0.02%, 0.07%, 0.01% for median values, as estimated from 
(X 2 - (X) 2 ) 1/2 /VN with N = 2 x 10 4 . 


and 20 keV. Ref. [T[ recommends 5(0) = 20.8 ± 1.6 eV b 
(quadrature sum of theory and experimental uncertain¬ 
ties). Our 5(0) is consistent with this, but the uncer¬ 
tainty is more than a factor of two smaller. Ref. [T| also 
provides effective values of S'/S = —1.5±0.1 MeV -1 and 
S"/S = 11 ± 4 MeV -2 . These are not literal derivatives 
but results of quadratic fits to several plausible models 
over 0 < E < 50 keV, useful for applications. Our values 
are consistent, considering the large higher derivatives 
(rapidly changing 5") left out of quadratic fits. 

The important quantity for astrophysics is in fact not 
S(E) but the thermal reaction rate; derivatives of S(E) 
are used mainly in a customary approximation to the 
rate integral mmm- By using our S' and 5" in a Tay¬ 
lor series for S(E) about 20 keV, then regrouping terms 
and applying the approximation formula, we find a rate 
(given numerically in the supplemental material) that 
differs from numerical integration of our median S(E) 
by only 0.01% at temperature T 9 = T/(10 9 K) = 0.016 
(characteristic of the Sun), and 1% at Tg = 0.1 (relevant 
for novae). 

How accurate is NLO? —Our improved precision for 
5(0) is achieved because, by appropriate choices of its 
nine parameters, NLO Halo EFT can represent all the 
models whose disagreement constitutes the 1.4 eV b un¬ 
certainty quoted in Ref. [lj- including the microscopic 
calculation of Ref. [5]. Halo EFT matches their S(E) 
and phase-shift curves with a precision of 1% or better 
for E < 0.5 MeV, and thus spans the space of models of 
El capture in the LER ITS] . 

The LO curve shown in Fig. [3] employs values of Ci, 
C 2 , ai, and 02 from the NLO fit. It differs from the 
NLO curve by < 2% at E = 0, and by < 10% at E = 
0.5 MeV. This rapid convergence suggests that the naive 
estimate of N2LO effects in the amplitude, (fc/A) 2 ss 4%, 
is conservative. And indeed, we added a term with this 
fc-dependence to the model, allowing a natural coefficient 
that was then marginalized over, and found that it shifted 
the median and error bars from the NLO result by at 
most 0.2% in the LER. Finally, we estimate that direct 
E 2 and Ml contributions to 5 in the LER are less than 
0.01%, and radiative corrections are around 0.01%. 

Summary — We used Halo EFT at next-to-leading or- 
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der to determine precisely the 7 Be(p, 7) 8 B S-factor at 
solar energies. Halo EFT connects all low-energy mod¬ 
els by a family of continuous parameters, and marginal¬ 
ization over those parameters represents marginalization 
over all reasonable models of low-energy physics. Many 
of the individual EFT parameters are poorly determined 
by existing S-factor data, at E > 0.1 MeV, but these data 
constrain their combinations sufficiently that the extrap¬ 
olated S'(20 keV) is determined to 3%. We estimate that 
the impact of neglected higher-order terms in the EFT 
on S'(O) is an order of magnitude smaller than this. 

Extension of the EFT to higher order and inclusion of 
couplings between s- and d-wave scattering states is not 
expected to reduce the uncertainty, although it would 
provide slightly greater generality in matching possible 
reaction mechanisms. There is, however, no indication in 
the literature that coupling to d-waves is important for 
S'(E) [5] in the LER. Our analysis could perhaps be ex¬ 
tended to higher energies, but for E > 0.5 MeV, accurate 
representation of Ml resonances is at least as important 
as reliable calculations of the El transition. 

The most significant source of uncertainty in our ex- 
trapolant is, in fact, the 1 keV uncertainty in the 8 B 
proton-separation energy, which can shift S (20 keV) by 
approximately 0.75%. This could be eliminated by better 
mass measurements. Further significant improvement in 
f>(20 keV) for 7 Be(p,y) 8 B requires stronger constraints 
on EFT parameters. Better determinations of s-wave 
scattering parameters seem to be of limited utility. The 
ANCs affect the very-low-energy S'-factor the most, and 
so more information on them, from either ab initio theory 
or capture/transfer data, would be useful. 

A number of other radiative capture processes whose 
physics parallels 7 Be(p, 7) 8 B are important in astro¬ 
physics. The formalism developed herein should be ap¬ 
plicable to many of them. 
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0 

'Z™ 

? 

0(3 Si) ( fm ) 

03 Sl ) (fm) 

ei 

Li (fm) 

C 2 5P9) (fm- 1 ) 

0(5 s 2 ) (fm) 

H 5 s 2 ) ( fm ) 

L 2 (fm) 

0.2336 

24.44 

3.774 

-0.04022 

1.641 

0.3269 
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SUPPLEMENTAL MATERIAL 


Common-mode errors for experimental data 


The quoted common-mode errors for Junghans et al ., sets BE1 and BE3, [37], Filippone et.al., [38], Baby et.al, 
[39l 40], and the Hammache et.al., 1998 data set [41] are 2.7%, 2.3%, 11.25%, 5%, and 2.2%, respectively. The 
data of Ref. [J2] are a measurement of the absolute ,15(186 keV) and of the ratios 5(135 keV)/5(186 keV) and 
5(112 keV)/5(186 keV). We treat each of these three quantities as one data point, so they do not need a CME. 


EFT details 

The modified-effective-range expansion for s-wave 7 Be-p scattering is: 

p(cot6 s (E) -i) 1 = + ^r s p 2 -2k c H{rj). (3) 

Here H(ij) = il>(ir]) + l/(2irj) — ln(i^), rj = kc/p , kc = ZrQ e Z p a em mu with mu the 'Be-p reduced mass, p = y/2mu.E, 
and 0 the digamma function m- 

An EFT parameter set that gives a good fit—as mentioned in the main text—is listed in Tabic [TT] 










7 


E (MeV) 

Median (eV b) 

—a (eV b) 

+ct (eV b) 

0 . 

21.33 

0.69 

0.66 

0.01 

20.97 

0.65 

0.63 

0.02 

20.67 

0.63 

0.60 

0.03 

20.42 

0.60 

0.58 

0.04 

20.20 

0.57 

0.55 

0.05 

20.02 

0.55 

0.53 

0.1 

19.46 

0.45 

0.44 

0.2 

19.27 

0.34 

0.34 

0.3 

19.65 

0.32 

0.30 

0.4 

20.32 

0.35 

0.31 

0.5 

21.16 

0.42 

0.41 


TABLE III: The median values and 68% interval bounds for S in the energy range from 0 to 0.5 MeV. At each energy point, the 
histogram of S is drawn from the Monte-Carlo simulated ensemble and then is used to compute the median and the bounds. 


Results for S'-factor and thermal reaction rate 

The median values and 68% interval bounds for S in 10 keV intervals to 50 keV and then in 100 keV intervals to 
500 keV is listed in Tabic HID 

Regrouping the Taylor series for S(E ) about 20 keV into a quadratic and applying the approximations of Refs. [21150 
yields 


N A {av) = - 


2.7648 x 10 5 


Tr 


■2/3 


exp 


'-10.26256 \ 

, Tl' 3 ) 


x(l + 0.0406X 9 1/3 


0.5099Tg /3 —0.1449T g + 0.9397Tg /3 -l-0.6791Tg /3 ), (4) 


in units of cm 3 s -1 mol -1 , where N A is Avogadro’s number. Up to Xg = 0.6, the lower and upper limits of the 68% 
interval for S{E) produce a numerically integrated rate that is 0.969(1 + 0.0576Tg — 0.0593X|) and 1.030(1 — 0.05T g + 
0.0511X|) times that of Eq. ©■ At Xg > 0.7 energies beyond the LER, and hence resonances, come into play and so 
these results no longer hold. We know of no astrophysical environment with such high Xg where 'Be(p, t) 8 B matters. 









